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ABSTRACT 

This paper presents a novel method for identifying indenta¬ 
tions on the boundary of solid 2D shape. It uses the signed 
curvature at a set of points along the boundary to identify in¬ 
dentations and provides one parameter for tuning the selection 
mechanism for discriminating indentations from other bound¬ 
ary irregularities. An efficient implementation is described 
based on the Fourier transform for calculating curvature from 
a sequence of points obtained from the boundary of a binary 
blob. 

Index Terms — blob analysis, shape indentations, curva¬ 
ture, Fourier transform 

1. INTRODUCTION 

Blob detection is a common unit in image processing pipelines 
in which an image is reduced to a set of distinct solid con¬ 
nected components that are then analyzed further. Most 
image processing packages provide a number of useful ge¬ 
ometric and topological quantities that can be computed for 
blobs, such as their centroid, area, eccentricity, perimeter, and 
so on. In some applications it is important to look at proper¬ 
ties of the boundary of the blobs. In particular, indentations 
and irregularities in the shape boundary that may have spe¬ 
cific meaning in a given domain (e.g., irregularly shaped cells 
in microscopy, defective widgets in manufacturing, etc...). 
How does one go about characterizing these indentations and 
irregularities to count them, locate them, and measure them? 

In this paper a method is presented to analyze blob bound¬ 
aries to identify indentations based on deriving a parameter¬ 
ized curve enclosing the blob and using measures of curvature 
on this curve to determine where indentations are present. A 
trigonometric polynomial representation of the boundary 
curve is obtained via the Fourier transform, and the necessary 
derivatives to compute the signed curvature are obtained us¬ 
ing convenient properties of the Fourier transform that make 
such computations very simple to state and implement. This 
method requires one parameter, severity , for tuning what 
constitutes an indentation versus a simple irregularity in the 
shape. This parameter has a geometric interpretation related 


to the radius of curvature along the boundary giving inden¬ 
tation selection a natural length scale. The method will be 
briefly compared with other approaches to this problem using 
common algorithms such as the convex hull of the shape to 
demonstrate where this method is preferable. 

2. APPROACH 

The key tool used in this method is the interpretation of the 
boundary of the shape as a parametric curve from which 
one can compute the approximate curvature at a sequence of 
points. When the sign of the curvature changes (i.e., when the 
normal to the curve switches from pointing into the shape to 
outside, or vice versa), one can infer that an inflection point 
has occurred corresponding to an indentation beginning or 
ending. The steps of the algorithm are: 

• Derivation of a trigonometric polynomial representa¬ 
tion of the discrete blob boundary point sequence as a 
continuous curve parameterized by arc length from an 
arbitrary boundary starting position. 

• Optional low-pass filtering of the parameterized bound¬ 
ary curve. This step eliminates high-frequency oscilla¬ 
tions in the curve due to the discretization of the con¬ 
tinuous shape when pixelized. 

• Calculation of the first and second derivative along the 
curve by element-wise multiplication in the frequency 
domain. Reconstruction of the first and second deriva¬ 
tive in the spatial domain necessary to calculate the cur¬ 
vature along the boundary is obtained via the inverse 
Fourier transform. 

• Identification of inflection points along the curve at 
points of curvature sign change. The regions delimited 
by these inflection points are then filtered by radius 
of curvature to eliminate detection of non-indentation 
features. 

A prerequisite for this algorithm is an existing segmen¬ 
tation pipeline that yields binary blobs from images, as well 
as a method for extracting an ordered point sequence along 
the boundary of these blobs. Such methods exist in many 



off-the-shelf software packages (such as the MATLAB image 
processing toolbox). 

2.1. Related approaches 

Shape analysis is a widespread activity in image processing. 
The majority of techniques found in the literature or pub¬ 
lished software packages rely on working directly with the 
pixel boundary of a shape or a polygon approximation for ei¬ 
ther the shape or its convex hull. The use of the convex hull 
for providing an approximation for a shape from which other 
properties can be easily computed is popular in large part due 
to the algorithmic efficiency and widespread availability of 
convex hull algorithms in popular packages using the Quick- 
Hull method [HO. As discussed later in Section 15.21 these 
methods make significant assumptions about what constitutes 
a shape indentation limiting their utility. 

Approaches that do not make assumptions about how in¬ 
dentations relate to the convex hull of a shape are more ro¬ 
bust to diverse shapes. Alpha shapes Eiia and Delaunay 
tri angulations GO can be used to calculate a boundary of an 
image from a point set (e.g., a set of pixels identified along 
the boundary of a binary blob). These methods unfortunately 
still require heuristic tests to determine when an indentation 
occurs given that the decision as to whether or not an indenta¬ 
tion is present requires more information than can be derived 
from a single boundary line segment or vertex. 

The calculation of the curvature using a Fourier series to 
represent the boundary curve naturally takes into account all 
points along the boundary. As such, the calculation of cur¬ 
vature at any point along the curve is intrinsically informed 
by the properties of the entire curve removing the need for 
any local heuristics near the point where curvature is calcu¬ 
lated. This simplifies the algorithm by removing heuristics 
that would induce additional parameters and potential points 
of brittleness. 


3. CURVATURE CALCULATION 

Given a blob with a corresponding boundary pixel sequence 
B = ..., ( x n ,y n )), the goal is to calculate the cur¬ 

vature k at each point along B. The signed curvature can be 
defined at the ith point along the curve in terms of the local 
first and second derivatives at that point: 




Fig. 1. Illustration of a parameterized curve boundary and the 
corresponding curvature measure along the boundary (clock¬ 
wise, starting at indicated point on boundary.) 


function represented as a Fourier series. Given some function 
x(t), the following property holds for its derivatives: 
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Where X(iuj ) = .F[a;(i)]. This allows the function to be 
differentiated in the frequency domain via simple multiplica¬ 
tion. In the case of blob boundaries, the function x(t) can 
be replaced by the sampled boundary points x\t]. The nth 
derivative in the spatial domain necessary for the calculation 
of k is recovered via the inverse Fourier transform: 


x {n) (t) = T~ x {(iuj) n X{ioj)) (3) 

This approach is attractive for a number of reasons. First, 
the Fourier representation does not use a spatially local ap¬ 
proximation for the curve in calculating the derivatives. This 
eliminates the need to select parameters such as a stencil size 
when using finite difference approximations for calculating 
derivatives directly from the boundary point sequence. Sec¬ 
ond, once the boundary curve has been moved to the fre¬ 
quency domain additional processing is possible to address 
issues such as high-frequency oscillations along a boundary 
due to pixel-level effects. Application of low-pass filters is 
very simple in the frequency domain and composes cleanly 
with the calculation of the derivative. In fact, any additional 
processing of the shape boundary that can be represented via 
convolutions can be included in the calculations applied to the 
boundary while in the frequency domain. 
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The derivatives along B must be calculated from the point 
sequence that was derived from the blob boundary. A simple 
and fast method for achieving this is via the Fourier trans¬ 
form (6|. The Fourier transform has a very convenient prop¬ 
erty that allows for simple calculations of derivatives of any 


4. INDENTATION IDENTIFICATION 

Identifying indentations on a blob boundary using curvature 
requires searching for inflection points where the normal to 
the curve changes sign. This occurs at points where the cur¬ 
vature reaches zero. Consider the example blob in Figure Q] 
The curvature along the curve is shown in the right-hand plot, 
showing clearly the points along the curve where the curva¬ 
ture passes through zero. 








Fig. 2. An illustration for the radius of curvature for different 
sized indentations. 

One parameter is used to tune indentation selection, sever¬ 
ity (a). This parameter has a close relationship with the radius 
of curvature that is used to define what constitutes an inden¬ 
tation. The magnitude of curvature, |re|, at a point along the 
boundary can be related to the radius p of the circular arc that 
fits the curve at that point 0: 

R“'’ (4) 

In regions of very low magnitude curvature, the corre¬ 
sponding circular arc would have a very large radius. Sim¬ 
ilarly, when the magnitude of curvature is high, the best-fit 
circular arc would have a low radius. Indentations can be 
thought of as places on the boundary where the sign of cur¬ 
vature is negative and the magnitude of curvature is high, and 
the use of the radius of curvature can provide a cutoff length 
scale for filtering indentations from non-indentation irregular¬ 
ities in which the sign of curvature is also negative. 

Given a maximum radius of curvature for what would be 
considered an indentation, p a , then p < p a can be interpreted 
with respect to curvature by Eq. Q]as |ft| -1 < cr -1 . Thus 
indentations will be considered only when |k| > cr. 

5. ALGORITHM 

The algorithm for the curvature calculation can be specified 
in MATLAB as follows: 

function k = curvature(xs, ys) 
fxs = fft(xs); fys = fft(ys); 

nx = length(xs); 

hx = ceil(nx/2)-1; 

ftdiff = (2i*pi/nx)* (0:hx); 

ftdiff(nx:-1:nx-hx+1) = -ftdiff(2:hx+1); 

ftddiff = (-(2i*pi/nx)"2)*(0:hx); 

ftddiff(nx:-1:nx-hx+1) = ftddiff(2:hx+1); 

dx = real(ifft(fxs.*ftdiff')); 
dy = real(ifft(fys.*ftdiff')); 
ddx = real(ifft(fxs.*ftddiff')); 
ddy = real (ifft(fys.*ftddiff')); 

k = (ddy.*dx-ddx.*dy) ./ ((dx."2+dy."2) . A (3/2)); 
end 


First, the x and y coordinate arrays of the points along the 
curve are mapped to the frequency domain with the FFT. The 
differential operator is calculated based on the length of the 
point sequences based on Eq. [2] This operator is then applied 
by pointwise multiplication against the frequency representa¬ 
tion of the points and the first and second spatial derivatives 
are then obtained by taking the real component of the inverse 
Fourier transform of the products (Fq. [3}. The calculation of 
k then follows exactly as defined in Eq.Q] 

Identification and counting of indentations then requires 
processing of the curvature sequence, as performed in the fol¬ 
lowing pseudocode. 

function n = countindentations(k, sigma) 
ksigns = sign(k); 
ksignchanges = diff(ksigns); 
idxchanges = [1; find(ksignchanges) ; ... 

length(ksignchanges) ] ; 

n = 0; 

for i = 1:length(idxchanges)-1 

ks = k(idxchanges(i):idxchanges(i+1)); 
if mean(ks) < 0 && max(abs(ks)) > sigma 
n = n + 1; 
end 
end 

return n 
end 

This algorithm takes the curvature sequence k and iden¬ 
tifies where sign changes occur. This is used to partition the 
curvature sequence into a set of regions with the same sign 
curvature. Each is then examined to identify those that repre¬ 
sent negative curvature regions in which the maximum curva¬ 
ture exceeds the value of the severity parameter sigma. 

5.1. Demonstration 

The application of this algorithm is shown in Figure [3] for six 
different shapes: two synthetic blobs drawn by hand, and four 
blobs derived from photographs of different types of flowers. 
In each of the sub-images the regions of negative curvature 
are indicated by the heavy boundary curve, and the title indi¬ 
cates the number of indentations that were detected. In some 
cases the effect of the cr parameter is apparent where irregu¬ 
larities along the surface were ruled out as indentations since 
their maximum curvature did not exceed the threshold cor¬ 
responding to the acceptable range of radii of curvature for 
indentations. 

5.2. Comparison with convex hull approaches 

A simple approach to solving this same problem using read¬ 
ily available algorithms is to examine the convex hull of the 
shape and count the number of regions where a gap appears 
between the convex hull and the curve itself. An example is 
shown in Figure[4]a), with the blob and its convex hull shown 




Fig. 3. Demonstration of the algorithm for six test images: 
two synthetic and four obtained from segmented photographs. 


in Figure |4}b). Gaps between the shape and the convex hull 
can be obtained by subtracting the blob from its convex hull as 
shown in Figure @tc). The problem with this approach is that 
it makes a simplistic assumption about what constitutes an in¬ 
dentation: specifically, that a shape without indentations will 
strictly match its own convex hull. This is unrealistic in many 
applications in which non-indentation boundary irregularities 
are expected. A consequence of this is that some indentations 
may be missed: a detailed region along the shape boundary 
is shown in Figure[4jd). Using the calculated curvature along 
the shape boundary is insensitive to this kind of shape since 
the curvature analysis makes no assumptions about the shape 
relative to a canonical analogue such as the convex hull. 

The attraction of convex hull methods is the ready avail¬ 
ability of efficient algorithm implementations. Fortunately, 
the method presented in this paper also leverages a fundamen¬ 
tal algorithm with widespread deployment, the Fast Fourier 
Transform. As such, the methods described in this paper are 
competitive in both accessibility as well as algorithmic ef¬ 
ficiency since the FFT has time complexity (){n log n) for 
n points, comparable to the average time complexity of the 
QuickHull method (0(nlogn) average, 0{n 2 ) worst-case). 

6. CONCLUSIONS 



(c) (d) 


Fig. 4. An example showing indentations missed by examin¬ 
ing the difference between the shape and its convex hull. 


via the Fourier transform and selection of indentations using a 
single scalar parameter corresponding to the maximum radius 
of curvature allowed to be considered an indentation. Unlike 
methods based on the convex hull or piecewise linear poly¬ 
gon approximations, the selection of indentations requires no 
heuristic tests and cleanly integrates frequency-domain filter¬ 
ing of the shape boundary to eliminate pixel-level discretiza¬ 
tion effects along the boundary of the shape. 

An implementation of this algorithm is available as an 
open source MATLAB package]]]. 
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